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Abstract. A class of diffusion driven Free Boundary Problems is considered which is character- 
ized by the initial onset of a phase and by an explicit kinematic condition for the evolution of the free 
boundary. By a domain fixing change of variables it naturally leads to coupled systems comprised of 
a singular parabolic initial boundary value problem and a Hamilton-Jacobi equation. Even though 
the one dimensional case has been thoroughly investigated, results as basic as well-posedness and reg- 
ularity have so far not been obtained for its higher dimensional counterpart. In this paper a recently 
developed regularity theory for abstract singular parabolic Cauchy problems is utilized to obtain the 
first well-posedness results for the Free Boundary Problems under consideration. The derivation of 
elliptic regularity results for the underlying static singular problems will play an important role. 
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1. Introduction. In this paper we consider a class of one phase Free Boundary 
Problems (FBP) characterized by the initial onset of a phase. Such Free Boundary 
Problems arise in the description of diffusion in polymers for instance. Under phys- 
ically natural conditions, these problems lead to a formulation in which the phase 
is initially absent. This feature manifests itself mathematically in that some of the 
equations in the nonlinear system become singular (if they are written in a fixed ref- 
erence domain) . In spite of the fact that this type of problems have been intensively 
studied in the literature over an extended period of time, all but one publication do 
not deal with the singular case considered here in more than one space dimension. 
The simplifying assumption that the phase be initially non-empty is typically added 
to avoid the mathematical complications stemming from the singularity. The one 
dimensional case has, however, been thoroughly investigated ([?, ?, ?, ?, ?, ?]) in 
the specific context of diffusion in polymers and, more generally, for diffusion driven 
FBPs ([?],[?] for instance). The methods used for the one dimensional case rely on 
the explicit use of the heat kernel to reduce the problem to the boundary. This ap- 
proach cannot be used in higher dimensions because the singular behavior induced 
by the initial condition on the FBP can not be decoupled from the diffusion oper- 
ator. This is due the the fact that the free boundary has a non trivial geometry 
in this case and to the fact that fundamental solutions (evolution operators) have 
not been studied for singular parabolic problems for which the singularity affects the 
underlying elliptic operator in an anisotropic way. The higher dimensional problem 
has recently been studied in [?] but only in its simpler quasi-stationary form. [?] 
established well-posedness for the problem in a class of functions which needs to be 
carefully crafted and leads to an asymptotic expansion for the singularity which is 
valid in the corresponding topologies and has practical relevance (cf. [?, ?]). Regu- 
larity results, maximal regularity in particular, for singular parabolic equations play 
a crucial role in this paper and have been recently obtained by the author in [?] 
where a construction of the evolution operator is given for a wide class of singular 
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parabolic problems. Maximal regularity is needed because the full problem will be 
solved by reduction to a model problem and perturbation about it. Since the ensuing 
perturbation is of maximal order (both in the singularity and differentiation) optimal 
regularity results are necessary. Previous abstract results were obtained by [?, ?], 
but they do not apply to the situation considered here because their validity does 
not cover the case of spaces of classical point-wise regularity, nor the case in which 
the singularity affects the equations anisotropically. Classical point-wise regularity is 
needed in the analysis proposed here since the singular parabolic equation is coupled 
to a Hamilton- Jacobi type equation modeling the front's dynamics. It has long been 
observed, as mentioned above, that one dimensional techniques can not be extended 
to the higher dimensional case. Summarizing, the presence of singular coefficients and 
the coupling to a Hamilton-Jacobi equations are two of the characterizing features of 
the problem under consideration. They make its analysis more difficult and delicate 
than that of the related but different classical Stefan-problem. This paper therefore 
offers a successful approach that fully overcomes these difficulties. 
The unknowns of the problem are a function u : Clt —>-M. defined on the open domain 
fit C ]R"+^ 3 {x,y) and its unknown boundary F*. The function u measures the 
concentration of the penetrant in the case of diffusion in polymers. Many geometries 
can be chosen for the domain Qf Here a strip-like setting is chosen where the domain 
is bounded by a fixed lower, and an upper moving hyper-surface denoted by Fq and 
Ft, respectively. Other configurations, like annulus-type domains are possible and 
interesting and the results obtained here would apply to those since they all would 
lead to the same local model problems. The system of equations satisfied by (u,Ff) 
proposed in [?] and generalizing [?, ?] then reads 



eut-Ax,yU = 0, in(J{t}xl]t, (1.1) 

u = g, on (0,oo) X Fo, (1.2) 



-d^,u = {I + eu)V , on[j{t}xTt, (1.3) 

t>o 

V = {l + 6Ut)u, on|J{t}xFt, (1.4) 

t>o 

f4^q=Fo, ati = 0. (1.5) 



In the above equations V denotes the front speed in normal (outward) direction Vt, 
and Ht its mean curvature. The function g > describes the concentration pro- 
ffie of a reservoir of the diffusing molecule whereas the positive constants e,d arise 
in the nondimensionalization process. They measure the deviation from a threshold 
concentration of the penetrant and the strength of the curvature term, respectively. 
Condition p. 3(1 models conservation of mass across the free boundary, whereas con- 
dition (|1.4|) is a phcnomcnological law which is capable of capturing the behavior 
observed in experiments know as case II or anomalous diffusion. The curvature term 
naturally appears if one assumes the front speed to be proportional to the local aver- 
age concentration in a reference ball of radius S > rather than to its point- wise value. 
The system is clearly nonlinear and, even if the boundary conditions were linear, it 
still would be. Indeed, two different solutions have different domains of definitions 
and cannot be added. This typical feature of FBPs becomes apparent after the simple 
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change of variables 

y 

{t, v) = {t, X, J , u{t, ^, 1]) = u{t, ^, ?7s(r, £,)) 

which transforms the problem to a corresponding one on a fixed domain. To do so, it 
is assumed that the fixed boundary of the unknown domain is given by 

To = R""^ X {0} 

and that the unknown moving boundary can be described as the graph of a func- 
tion s(t, •). Latter assumption is motivated by (initial) condition p.5|l . In the new 
variables, the domain becomes the strip 

5 = M"-i X [0,1] 

with the obvious boundaries denoted by Fj, j = 0, 1. Using the old notation for the 
new variables the system now reads 

l + j/2|Vs|2 5 1/1.^ N sAs-2|Vs|2 
£Ut - L\xU 5 o„u = ey-OyU ~ 2y-[Vs\OyVu) - y dyU 

(1.6) 

u = g on To (1.7) 

1 + |Vs|2 

— -dyu={l + eu)s^ {Vs\Vu) onTi (1.8) 

s = Vl + |Vs|2(l + ,5Hf)u onTi (1.9) 
s(0, •) = (1.10) 

In [?] the quasi-stationary approximation (e = 0) in two space dimensions was consid- 
ered in the presence or absence of the mean curvature term in two space dimensions. 
Here there is no restriction on the spatial dimension and the full evolutionary problem 
is analyzed but the curvature effects are neglected (6 = 0). 

The main result of this paper is establishing the well-posedness of system (|1.6ll - (|1.10|l . 

To do so, appropriate function spaces have to be used which are able to capture the 
regularity of the solution as well as its asymptotic behavior at the origin. The choice 
of function spaces is limited by the simultaneous presence of a singular parabolic 
problem and a Hamilton- Jacobi equation. Spaces of classical regularity are better 
suited for the latter and thus reduce the freedom of choice for the parabolic problem. 
A compromise can be reached by using spaces of classical Holder regularity in space 
and of singular Holder behavior in time. A number of preparatory results are needed; 
they are formulated in the next section. In section 3 elliptic and parabolic results 
for singular equations will be derived which play a crucial role in the existence proof 
given in the last section. 

2. Preliminaries and Setting. In order to obtain existence results for system 
H1.6|l - (|1.10|l a linearization procedure will be used in combination with maximal reg- 
ularity results for the relevant linearized problems. The presence of the singularity 
complicates the analysis significantly since the relevant regularity results are not avail- 
able and need to be derived first. The linearization procedure is meant to capture 
the leading order terms of the differential operators at the origin both with respect 
to their differentiation order and their singularity degree. 
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Of interest in this paper are classical solutions. If such a solution existed, one would 
be able to guess from (fr7|) and lfO|) - ((T~Tn|) that, 

s{t, x) w tg{x) for t ^ . 

This indicates that the relevant model problem which captures the leading order 
behavior of (|1.6|I - H1.8|I has the form 

eut - A.^u - -jj^^j^dlu = f{t, X, y) in (0, oo) x , (2.1) 
u^g{t,x) on(0, oo)xro, (2.2) 
jdyu = h{t,x) on(0, oo)xri, (2.3) 

where in the situation at hand g is independent of the time variable and c — g. The 
more general case of c independent of g is, however, of interest and motivates the 
different notation. For reasons which can be guessed and will become apparent later, 
both parabolic (e = 1) and elliptic (e = 0) regularity results are needed in the analysis 
of l|2.1|l - (|2.3|l . The main ingredients needed to derive such regularity results are 
vector-valued Fourier multiplier theorems and the use of spaces of singularly Holder 
continuous functions in the time variable and the construction of an evolution operator 
for singular families of generators. The first are crucial in the analysis of the elliptic 
problem with time frozen and the latter allow for the quantitative characterization of 
the singular behavior in the origin {t = 0) both in the elliptic and the parabolic case 
(where maximal regularity results are needed). 

For the sake of completeness the formulation of the relevant Fourier multiplier theorem 
and the definition of the classes of singular Holder continuous functions needed in 
the analysis are given here. The basic observation illuminating the reason for their 
combined use will also be presented in this section. 

Assume that E is a. given Banach space and that T > and /3 e (0, 1) [U{1 — }], then 
the standard Holder [Lipschitz] space is given by 

C^{[0,T],E) {/£ C([0,r],i?) I M^,[oT] := sup M^L.^ < ^} 
with norm 

II ■ 11/3 = II • l|oo,[0,T] + [-l/SJCT] ■ 

If /? = 1 is chosen, one obtains the space of Lipschitz continuous functions. To 
distinguish it from the space of continuously differentiable functions the notational 
device /3 = 1— is used. This means that (3=1 for all practical purposes except in the 
notation for the space which becomes C^^. Singular counterparts are given by 

C^((0,r],i?) := {/ e B{{Q,T],E) \ [t ^ t^u{t)] e &{{Q,T],E)} (2.4) 

with weighted norm defined through 

II • h,l3 ||w||oo,(0,T] + [(•)^u]/3,(o,T] • 

The symbol denoting the time interval in the notation for the norm will be dropped 
in the sequel with the understanding that the interval of definition does not contain 
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the origin for singular spaces, whereas it does for regular ones. The following closed 
subspace of regular Holder functions will also be useful 

C^([0,T],£;) := {/ e C^{[0,T],E)\f{0) = O} . 

The vector- valued Fourier multiplier theorem which is needed here can be found in [?, 
Theorem 6.2]. We reproduce it here for the ease of the reader. The following symbol 
class 

S"'{Eo,Ei) {a e 0"+^ (K" \ {0}, /:(i?o, | 

sup |[l + |en('"+l"l)/29M0|^(s„^.<oo, \a\<n+l} (2.5) 

where Eq and Ei are Banach spaces and m S Z, is instrumental in the formulation of 
the result. Given a symbol a € S"^{Eo, Ei), an operator can be associated to it by 

a{D) := T^^aT 

through conjugation with the vector-valued Fourier transform 

T e A. S,)) n (5'(M", S,)) , J = 0, 1 , 

where S and S' are the operator- valued Schwartz' spaces of fast decaying test functions 
and of tempered distributions, respectively. They are endowed with their natural 
topology. 

Theorem 2.1. /?, Theorem 6.2] Suppose that B G {B,B,h] and that m £ M. 
Then 

[a ^ a{D)\ e c{s"\Eo,E,),C{B;,,iW^, Eo),B;+^(W^, E,))) 
for s e M and p,q e [1 , oo] . 

o 

The symbol B denotes the regular Besov spaces whereas B^ ^ and bp g denote the 
closures of S and B^'^^ in Bp ,^, respectively. Various equivalent definitions are given 
in [?]. Here only regular Besov spaces with p = q = oo of positive fractional order 
s G \ N are used in which case one has 

B^,^=BUC^(R",£;), 

where the space on the right is the standard space of bounded and uniformly Holder 
continuous functions given by 

BUC"(M",£;) {/eBUCt"l(M",£;)|a"/eBUC""t''l(M",£;)V|a| < [s] } . 

A crucial observation connecting the dilation of symbols of type H2.5|l and singular 
Holder spaces was already obtained in [?]. It allows one to deal with singular elliptic 
and, eventually, singular parabolic boundary value problems. 

Proposition 2.2. [?, Lemma 2.5] Assume that a E S"^{Eo, Ei) for some m G Z. 
Then 

[t^t^ata] eC}-((0,T],S™(£;o,i?i)) 
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for at(a) := a{t-). 

Remark 2.1. A straightforward adaptation of the proof given in [7] of Proposition 
also shows that 
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^rVta] e ci-([o,r],s™-i(i?o,i?i)) 

The importance of this proposition Ues in the fact that muhiphcation of singularly 
Holder continuous functions is well-defined and continuous as a map defined in various 
combinations of spaces 

(2.6) 
(2.7) 
(2.8) 
(2.9) 

for a,/3 e (0, 1) U {1 — } and (3 < a. The proof is elementary and can be found in [?] 
along with the (natural) definition of the spaces for negative lower indices. 
Finally optimal regularity results are yet another essential ingredient to deal with the 
full nonlinear problem. It provides maximal regularity results for the relevant class 
of singular parabolic problems in classes of singularly Holder continuous functions 
which are, as they need to be, perfectly compatible with the corresponding results for 
singular elliptic problems in the same class of functions. The result needed here has 
been derived in [?] and is formulated in the next theorem. It gives conditions for the 
well-posedness of the singular abstract Cauchy problem 

u- A{t)u^ f{t),t>Q (2.10) 

for a Banach space valued function u : (0,r] Eq and a singular family of "elliptic 
operators" A (that is, of generators of analytic semi-groups). The symbol H~{Eo,u}) 
denotes the class of generators of exponentially decaying semi-groups which are not 
necessarily strongly continuous (as introduced by [?]). 

Theorem 2.3. /?, Corollary 3.3] Assume that A satisfies the following assump- 
tions 

(i) A{t) eH-iEo,uj), t>0, (2.11) 

(n) \\[A{t)^A{s)]A-\r)\\c^Eo)<c^-^ and 

II [A{t) A{s)] i-Ar^ir) WciEo) < c{t - s) , (2.12) 
(iii) lim A-\t) =0, (2.13) 

for some pe (1,2) and < t < s < t < T . Let f e C'^{{0,T], Eq) for some f] e (0,1) 
and ^ = 0,/3. Then (|2.10|) has a unique solution u € C^^((0, T], i?o) satisfying 

u,Au e Cf^ and ||u||/5,^ + ||Au||/3,-y < c||/||/3^^ 

where \\ ■ ||^,o = II ■ 11/3- 

Using a combination of elliptic and parabolic estimates it will be possible to obtain 
a satisfactory regularity theory for H2.1|l - (|2.3|) . The latter will used in the analysis of 
the full evolution system H1.6|) - (|1.10|) for 6^0. 
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3. Elliptic and Parabolic Estimates. In order to use the results presented in 
the previous section in the analysis of (|l.t)|l - H1.10|) . appropriate function spaces have 
to be chosen in which to work. The choice can typically be justified by balancing the 
regularities in the nonlinear equations in such a way that a fixed point argument can 
be applied to the set of equations. Whereas a variety of function spaces of Sobolev 
type are available for dealing with elliptic and parabolic problems, the fact that the 
system of interest contains a Hamilton- Jacobi equation for the evolution of s, which, in 
turn, appears in differentiated form as a coefficient in the equations for the evolution 
of u restricts the choice to spaces of classical regularity for s. It is indeed impossible 
to work in Sobolev-Slobodeckii spaces for u as the loss of regularity incurred in taking 
traces on the boundary cannot be made up for by H1.9|l which does not possess any 
regularizing effect. It becomes clear that spaces of classical regularity for both u and 
s need to be chosen. In order to use the Fourier multiplier Theorem 12.11 the choice 
is therefore reduced to spaces of bounded uniformly continuous functions in the x- 
variable. As far as the y variable is concerned a good choice is given by the space of 
continuous functions as will become clear later. The "base" space Eq for u is therefore 
chosen as 

£;o BUC^+"(M",C(0,1)) (3.1) 

for a £ (0, 1) and where C(0, 1) denotes the standard space of continuous functions on 
[0, 1]. The space with one less regularity degree in x would seem like a more natural 
choice but the choice made here makes it much easier to deal with inhomogeneous 
Neumann boundary conditions such as (|1.8() . This point will become clear later in the 
analysis in Section rOl In view of the singular non-autonomous nature of the prob- 
lem, the standard procedures ([?]) to deal with inhomogeneous (nonlinear) boundary 
conditions cannot be utilized. 

It follows that a proper choice of "base" space for the full evolutionary problem is 
given by 

Eo := C^((0,T],BUC^+"(R",C(0,1))) (3.2) 

for u and by 

§ := C^+''([0,T],BUC^+"(R")) n C''([0, T], BUC^+"(R")) (3.3) 

for s. A closer look at the underlying elliptic problem (e = 0) is necessary in order 
to better deal with the inhomogeneous boundary terms in H2.1|l - (|2.3|) . The operator 
— ^x — t'^Ji(x) '^^^ clearly non-constant coefficients and multiplier Theorem 12. II is 
formulated only in the translation invariant case. It therefore needs to be shown that 
localization arguments apply in the operator-valued setting considered here. 
The first step is, however, to obtain regularity results in the constant coefficient case. 

3.1. The Elliptic Case - Constant Coefficients. Consider problem (|2.1(l - 
H2.3|l and assume that £ = and that c is a constant which can be assumed to be 
1 without loss of generality. Then taking a Fourier transform in the a;- variables, one 
obtains a parameter dependent boundary value problem for u which can be solved 
explicitly. The solution has the structure 



u{t, e, ■) = iVtc(e, •)/ + t<^tb{t ■% + (Tta{^, -Yg (3.4) 
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where at denotes dilation by t > and the operator-valued symbols a, b, c are given 

by 

cosh(l^l) l^lcosh(l^l) 

The first two symbols have to be considered as multiplication operator valued with 
respect to the variable y, whereas the last contains the sectorial operator C which is 
the operator —dy on C(0, 1) with domain 

dom(C) = {ue C^(0, 1) I u(0) = , dyu{l) = 0} =: Co_o(0, 1) ■ (3.6) 

The spaces 

:= {u G C''(0, 1) I u{0) = 0} , /c = 0, 1 

will also be needed. Proposition 12 . 21 can be used in combination with the multiplica- 
tions H2.6|l - H2.9|l to obtain regularity results in weighted Holder spaces provided (|3.5|l 
are in the appropriate symbol class. 

Theorem 3.1. The symbols (|3.5|1 satisfy 

ae50(C(0,l),C(0,l)), 

6e5'=(C(0,l),Cj-'=(0,l)), 

ce5^"(C(0,l),Cg7o^(0,l)) 

for fc = 0, 1 and j = 0, 1, 2. The notations for the spaces have to be interpreted so that 
boundary conditions are imposed only if enough regularity is available but otherwise 
disregarded. 

Proof. Take the last symbol first. It is known that the operator C is sectorial and 
invertible (cf. [?]). It follows that 

ll(iei' + C)-i||£(c,c-^)< c^-j^, j = 0,l,2 (3.7) 

As for the derivatives, it is seen by induction that, for ^ 7^ 0, 

\0\<\a\ 

for polynomials Pa^fi of degree at most \f3\. The estimate therefore follows from H3.7|) . 
Next consider the "boundary symbols" a and b. A simple series expansion reveals 
that both symbols are actually analytic functions of (and y of course). It follows 
that they are smooth for all ^ G M" in spite of appearances. It is therefore legitimate 
to concentrate on the decay properties of the symbols without necessarily using the 
nowhere vanishing weight (1 + |CP)^^^ needed in the definition of the symbol class. 
Observe that both are smooth functions of y regardless of ^. Consequently they define 
multiplication operators for C^~'^(0, 1) and = 0, 1, 2. Also observe that, for b, the 
smoothing effect in x is due to the decay in ^. They clearly satisfy the required 
boundary conditions. They further satisfy the relations 

dybit y) = a(e, 1 - y) , dpi^, y) ^ v) , 

9,a(C, y) = -|ePfo(e, 1 - y) , dlaii, y) = |epa(C, y) , 
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which make it easier to derive their mapping and symbol properties from one another. 
For instance, if it were known that a, 6 G 5°(C,C) the first relation would imply 
that dyb e 5°(C, C) and therefore that b e 5°(C,C^). The last one, on the same 
assumptions, would show that a G 5~^(C,C^). More details are now given for the 
analysis of symbol a. Corresponding claims for the symbol b can be obtained by 
similar calculations. The claim is that a G <S*'(C, C). It is plain that 

sup sup \a{£_, y)! < oo . 

CeR" yG[0,l] 

As for the first derivatives in ^ one has 



dAi, y) = ^ ^^^^^(1^1) [(1 - y) sinh(|e|(l - y)) - tanh(l^l) cosh(|e|(l - y))] = 

|.a(C,2/)[(l - y) tanh(|^|(l - y)) - tanh(|e|)] . (3.8) 

The maximum in y is either taken on the boundary y = 0, 1 or in the interior. On 
the boundary it either vanishes or is exponentially decaying in |^|. It therefore only 
needs to be controlled in the interior to conclude its estimation. Only the term in 
brackets in the right-hand-side in the first line of H3.8|l depends on y. Setting its first 
y-derivative equal to zero, one arrives at the equation 

tanh(|e|(l - y)) _ 1 



(l-y)l^l |^|tanh(|e|)-l 
which, for |^| large, is asymptotic to 

1 1 

i^-ym^Wl' 

Thus the maximum is located at y ~ for large |^|. As for the maximal value for 
large |^| one has 



1 , ?j ,e 1 ^ 1 



5j-«(e,-)-^T^a(C,-)- (3.9) 
which entails the desired estimate 



sup sup {l + \^\^y/^\dja{^,y)\<^. 

xGK" i/G[0,l] 

Equation H3.9|) follows from 

(l-^)tanh(|C|-l)-tanh(|^|) (3.10) 

observing that tanh(s) = 1 up to exponentially small terms in s for s > large, that 
is, 

|tanh(s) - 1| < ce-^'s > 0. 
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The second derivative satisfies 



The first term in the right-hand-side can be seen to decay hke 1/|CP in view of the 
explicit pre- factor and the estimate obtained for the first derivative (|3.8|l . The second 
and third terms add up to 

|jf a(e, v) [(1 -vf-l- 2(1 - y) tanh(|e|) tanh(|C|(l - y)) + 2 tanh^d^D] 

Similar calculations as those performed for the first derivative show that the maximum 
is now attained at 

2 



and it amounts to 



2 . 2 2 4 



Latter follows from 



(1 - ll )' - 1 - 2(1 - li ) tanhdei - 2) tanh(|e|) -f 2 tanh^d^]) ^ A 

which uses (|3.10|) again. In conclusion it is obtained that 

sup sup (1 + |^n|9,9ja(^,?/)| < oo. 
5eR" iye[o,i] 

Comparing with the proof of [?, Lemma 2.2], it is observed that the arguments are 
almost identical. The only difference are the additional factors like for the first 

derivative or and for the second. Terms containing these can always be han- 
dled through their explicit dependence on ^ and the estimates from previous deriva- 
tives. The inductive argument used in [?] can therefore be adapted to the current 
situation and leads to the desired result. □ 

Combining these results with Theorem 12.11 Proposition 12 . 21 and H2.6|l - (|2.9(l the next 
important theorem is obtained. To simplify its formulation we introduce the abbre- 
viated notation 

C? BUC ■ C? 7 

for the function spaces 

C?((0,T],BUC-(R",C?,7(0,1))) 

where the question marks can be substituted by any of the relevant regularity and 
singularity parameters. 

Theorem 3.2. Let a,(3 € (0, 1) and — 0,(3. Given 



(/, g, h) e C^ BUC^+" C X Cf' BUC^+" x Cf' BUC 



2+a 
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there exists a unique solution u of (|2.1f) - ()2.3|) for c = const which belongs to the space 

d 

El ■- {u e Eo I d^{-^yu e Eo for = 0,1,2 s.t. \a\ + j < 2} 

where 



Eo := C^BUC^+"C . 

It is given by 

u = T^H^atcT f + T^^tatbT h + T^^atoT g . 

= :A(t)-l =:B.N{t) =:Rn{t) 

Modulo the singularity, the terms in the above representation clearly point to their 
asymptotic behavior at the origin. This result needs to be extended to the non- 
constant coefRcient case. This is done in the next section. 

Remark 3.1. To make it easier to refer back to the above result in the case 
7 = 0, the corresponding spaces will be denoted by E^.Q; fc = 0, 1. 

3.2. The Elliptic Non-constant CoefRcient Case. In this section we adapt 
an abstract formulation of a localization argument proposed by Angenent in [?] (see 
also [?, ?]) for M" to cover the case of boundary value problems in an operator- valued 
context. It is used in combination with the Fourier multiplier Theorem 12.11 to give 
the basic regularity results needed for the singular elliptic boundary value problem 
underlying H2.1(l - (|2.3|l in the case of constant coefficients. It should be observed that 
localization arguments rely on perturbation results for differential operators by lower 
order terms. These are small compared to the leading order operator only in a qual- 
itative sense. Perturbation results can therefore only be applied for the resolvents 
with large A where the decay properties of the inverse to the leading order operator 
yields smallness of the lower order terms. The structure of the singular operators 
considered here can be exploited in order to avoid shifting the leading order operator 
and in order to obtain direct invertibility results. The standard argument of course 
shows that non-constant coefficient operators are sectorial whenever their constant 
coefficient counterpart is (see [?]). This remains valid here and can be used if the 
time variable is kept fixed. 

The goal of this section is to extend the validity of Theorem 13.21 to the non- 
constant coefRcient singular Boundary Value Problem 

{A{t), 70, -^lidy) : El ^ Eo X ^Ei (3.11) 

where A(t) is the elliptic operator driving the evolution in (|2.1|1 and the "boundary 
space" 9Ei is defined by 

5Ei = BUC^+" X BUC^+" =: do^i x On^i , 7 = 0, . (3.12) 

Observe that Theorem 13.21 simply states that the above singular Boundary Value 
Problem has a bounded inverse in the given topologies and for constant coefficients. 
The localization procedure proposed in [?] is abstract and is based on the concept of 
resolution. A resolution of a Banach space E is simply a triple [F, e, 5) where F is 
a Banach space and the maps 6 : E F and e : F ~* E satisfy e o 5 = id^. The 
operator 5 plays the role of the localizing operator, whereas e resynthesizes the local 
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contributions. If the space E has subspaces of interest, they should correspond to 
"similar subspaces" of F and should be left invariant by the resolution maps. The 
localization operator can be constructed as follows 

6 -.Ej ^ ®kezEj , u ^ {uk)kez" iu(pk,r)kez^ , j = 0, 1 

where F*^ = Kj for k e Z" and the sequence space Qkez^j is endowed with the 
supremum norm sup^.^^ . The maps ((^^ r)k£Z^ chosen as to form a smooth 

resolution of the identity in M" 



E 



"Pk-, 



k£2 



subordinated to a cover by cubes Qk£Z",r of fixed side length r > obtained by 
translation and dilation of the standard cube of side size r = 1 centered at the origin. 
The support of ipk^r should contain Qk^r and be contained in the union of at most 
finitely many adjacent cubes. The synthesis operator is then given by 

£ : ©fcezF*" Ej , {uk)keZ" ^ "fc¥'fc,r 

in which case the relation defining a resolution is clearly satisfied. Observe that the 
mappings d and e are continuous. The idea is to define an operator on (Bkez^i which 
parallels the Boundary Value Problem (|3.11(l and to use its inverse to approximate 
the desired inverse of (|3.11(l via the use of the maps e and 6. Define the operator 
{A',B') on ©fcezFi through 

{A'',B'')uk = {-A^uk - ,2 ■. dyUk,joUk, ^jidyUk) 

where Cfc|supp(<^fc ^) = c\supp{ip^ ^) and Ck is otherwise smoothly extended without in- 
creasing its norm to the whole space (the details of the extension procedure can be 
found in [?, ?]). The Boundary Value Problem (A'^^B'') can be made arbitrarily close 
to a constant coefficient one, say by substituting Ck by Ck{xk) and making r small. It 
is therefore invertible and enjoys the properties claimed in Theorem l3.2l The diagonal 
operator {A',B') is then invertible itself and maps x 9Fi] to (Bkez^o where 

d¥i is defined in the obvious way. Denote by [A'{t)-^ R'jj{t) R'j^{t)] its inverse. 
Then 

s[A'{tr^ R'^it) R'^{t)]S (3.13) 

should represent an approximation to the solution operator for the singular boundary 
value problem H2.1|l - I|2.3|l . In order to show this, it needs to be shown that H3.13|l is 
an approximate left and right inverse. Starting with the latter, compute 



At)' 

70 



e[A'{t)-^ R'j^it) R'^{t)]5 

idv,+[A{t),e]A\ty^5 [A{t),e]R':j{t)5 [A{t) ,e]R'j^{t)5 

id^„Ei (3.14) 

id^ia^Ei 
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where [^(i), s] — A{t)e — eA{ty . Invertibility would follow if 

idE,+[A{t),e]A'{t)-^5 (3.15) 

were invertible. Here the specific structure of the operator A{t) needs to be exploited 
since standard perturbation arguments would only allow to show that the resolvent 
TZ{X, A) exists for large enough A, in which case the lower order commutator term is 
small. The commutator term (|3.15|) can be computed to give 

^ [2(V^fe,,|V-) + A^k,r-] - T^^TT^y] "Vfc/ ■ (3.16) 

To see that the inverse appearing in (|3.16|) does indeed exist on the desired spaces, 
observe that r > can be chosen such that 

is small in £(Ei , Eq). The realization A'' of A'^ with homogeneous boundary conditions 
is therefore invertible and the constant coefficient estimates of Theorem (|3.2|) carry 
over to it (with larger constants, of course). Now, the estimates 

|i2V^fe..V[-A, - -^l^dir'WciEo) < cT and 

show that the commutator is indeed small, at least for T small enough. Here the 
decay properties of the resolvent as t — > are used as a substitute for making A large 
in the standard argument. Heuristically this makes sense, since making t small, just 
as making A large, makes the operator "more and more elliptic" . This is due to the 
specific nature of the singular operator. Since the difficulty of the problem stems from 
the origin, making T small does not in any way weaken the result. 
Next it needs to be shown that H3.13() is also a good approximation for a left inverse. 
To that end, observe that 

u - s{A'it), B'{t))~^S{A{t),Bit))u 

is the same as 

u- (eA'{t)-^SAit)u + eR'D{t)6jou + eR'N{t)S-fi^uj ^ eA' {t)-^[Ait),S]u (3.17) 
because 

u = eA'{t)-^A'it)5u + eR'jj{t)5-iou + eR'j^{t)-fi^u 
by definition. The notation 

[Ait), 6] =A'{t)d~dAit) 
was used in (|3.17|l . It follows that 

e{A'{t),B'{t)y's{A{t),B{t)) - idg, -eA'{t)-^[A{t),S] 
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The term containing the commutator is lower order and can be estimated just Hke 
before exploiting the structure of A{t) in order to yield the invertibility of 

idE, -eA'{tr^[A{t),S] 

and, thus, a left inverse for the singular boundary value problem. The result is 
summarized in the next theorem. 

Theorem 3.3. Let a,/3 G (0,1), j = 0, P and < cq < c e BUC^+". Given 

(/, g, h) e BUC^+" C X BUC^+" x BUC2+" 

there exists a unique solution u of (I2.1|) - (|2.3|l which belongs to the space 

d 

El ;= {u eEo e Eo for |a|,j = 0,1,2 s.t. \a\ + j < 2} 

where 

Eo := C^BUC^+"C . 

It naturally splits into three components 

u = A{t)-^f + RN{t)h + RD{t)g 

with the same asymptotic behavior at the origin as in the constant coefficient case and 
where 

A{t)-'f = {A{t),B{t))-\f, 0,0) , RD{t)g = {A{t),B{t)y\0,g,0) 
and RN{t)h= {Ait),B{t)y^{0,0,h) 

respectively. 

Proof. The only part of the proof missing is for the claim about the asymptotic 
behavior of the solution. It is obtained using the representation 

"T -T[Ait),e]R'jj{t)6 ~T[Ait),e] R'j^{t)6 
id7oEi 

for the inverse of 

'idEa -[Ait), e]A'{t)-^S [As]R'jy{t)5 [A{t),e]R'j^{t)5 
id7oEi 

id-yia^Ei 

where T — [idE„ —[A{t),e]A'(t)^^S] ^, the factorization of the resolvent implied by 
and the mapping properties for the operators involved which follow from the 
symbol analysis combined with Theorem 12. II Proposition 12.21 and H2.6|l - H2.9|l . Take 
for instance the entry TA{t) e R'j2,{t)S. The claim follows from 

TeC{Eo), ^e/:(Ei,Eo) and eR'^S E £{80^1,^1) . 

□ 

From now on, the tildes on the solution operators of the non-constant coefficient 
case will be omitted. Since the corresponding operators in the constant coefficient 
case will no longer be used in the analysis, no confusion seems likely. 
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3.3. The Parabolic Problem. It is now possible to return to the analysis of 
the model problem H2.1|l - (|2.3|) for e = 1. If the boundary conditions are homogeneous, 
then Theorem l2 .31 gives existence and regularity of a solution u e Ei whenever / e Eg . 
The notation Ei (see Theorem 13.211 has been used so far to describe the regularity 
space of the singular elliptic boundary value problem. In the parabolic case the same 
notation indicates the space 

d 

El := {u G Eo I dtu e Eq , d'^{-^yu £ Eq for \a\,j = 0,1,2 s.t. |a| + j < 2} 
in accordance with Theorem 12 .31 To obtain optimal regularity results for 



it is therefore sufficient to check that conditions H2.11|l - (|2.13|l are satisfied. The op- 
erator C was defined in (|3.fci|) . Freezing coefficients arguments such as those in [?, ?] 
or [?] show that A{t) does indeed generate an analytic semi-group on for each 
fixed t > which is exponentially decaying since C is invertible (see also the begin- 
ning paragraph of Subsection 3.2). Estimates (|2.12() - (|2.13() follow from the regularity 
theory developed in the previous section and the fact that A{t) generates an expo- 
nentially decaying analytic semi-group. The latter makes it possible, in particular, to 
define the fractional power appearing in 1)2.12(1 (see [?] for more details). Condition 
(|2.13|l follows from 

II^W-'ll£(Eo) < CT^ 

which, in its turn, follows from the observation that (^^^jyC)^^ is bounded in Eg and 
from 

c^(a;) 

The first condition in (|2.12|) follows from 

[A{t)-A{s)]A{r)-^^ ^^-jy 4-\-r^A^ ' ^^"^ 



and the fact that 



C r 2. C 



c^(x) c^(x) 



£(Eo) 



is uniformly bounded for r G (0,r]. The second condition in H2.12|l follows similarly 
by using abstract mapping properties of fractional powers and, in particular, that 

It remains to be shown that problem 1(2.111 - 1(2.3(1 can be solved for inhomogeneous 
boundary conditions as well. The result is formulated in the next theorem. 
Theorem 3.4. Let c G BUC^+" and 

/ G Eo , g G BUC^+" n Cl+'^ BUC^+" , /i G B\JC^+°' n Ci+'^ BUG" . (3.18) 
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Then there exists a unique solution u £ Ei of (|2.1|l - (|2.3|l . 

Proof. It can be assumed without loss of generality that / = 0. Looking for a 
solution u in the form 

u = v + RD{t)g + RN{t)h 

one obtains that v satisfies 

V - A{t)v = ^ [RD{t)g + RN{t)h] . 
Using equations H2.1|) - (|2.3|) . it can be checked that 
d 2 

-[RD{t)g\ = -A{t)-^ ^^RD{t)g + RD{t)g 

^ [RN{t)h] - -^A{t)-^A^RN{t)h + -^RN{t)h + RN{t)h . 

The regularity results obtained in the previous section combined with the assumptions 
then imply that 

{^^Jt t^^^^)-*^ + RN(t)h] } e Eo . 

Take for instance \t ^ jA{t)^^l\xRD[i)g\ and observe that 

2 

Rd G £(9dEi,Ei) , e /:(Ei,Eo) and [t ^ e /:(Eo) 

or RD{t)g in which case the stated regularity follows from 

g e 7oEo and Rd £ >C(7oEo,Eo). 

The fact that Rd G 'C(7oEo,Eo) has not been explicitly proven but can be obtained 
by symbol analysis and freezing coefficients along the lines of Sections 3.1 and 3.2. 
The claim then follows from Theorem 12. 31 □ 

4. The Hamilton-Jacobi Equation. Consider the Hamilton- Jacobi equation 



St - Vl + I Vs|2 w = in (0, oo) X M" , (4.1) 
s(0, •) = on R" (4.2) 

for a given w G C'^ BUC^^" satisfying w(0, •) — g. The method of characteristics allows 
one to recast this Hamilton-Jacobi equation as a system of ODEs in the following 
manner 

i = l, i(0) = 0, (4.3) 

x = , ^ v(t,x), a;(0)=pGR", (4.4) 

• - w(t,x),z(0) = 0, (4.5) 



Vl + bp5tw(t,x), r(0)==.9(p) (4.6) 



P= Vl + H'V.w(t,a:), p(0)-0. (4.7) 
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This system is easily seen to reduce to 

X = j=— v{t, x) , x{0) = peR\ (4.8) 

P=yi + \p\^VMt,x),p{0)^0 (4.9) 

as all other unknowns can be obtained after solving this reduced system. The as- 
sumption on V makes it possible to solve this system on a possibly small time interval 
[0,T] which is independent of p G R". Exploiting the regularity assumption on v it 
can easily be seen that the flow mapping 

satisfies 

[t ^ {Xt,Pt)] e C''BUC^+"(M2",M2n>) 

where the map is obtained by solving system (|4.8I) - H4.9|I with the second initial con- 
dition substituted by p{0) = rj. It is clearly a flow of diffcomorphisms. Furthermore, 
since an equation satisfied by (DpX,Dpp) is easily derived from H4.8() - (|4.9() and since 

{Dpx{0, p, ri),Dpp{{), p, 7^)) = (idK„, 0) , 

it follows that DpXt {p, 0) satisfies 

IlidK. -DpXt{p,0)\\ < i 

uniformly in p G M". This implies that Xt{-, 0) is a diffeomorphism and yields uniform 
estimates for the inverse of Xt. It follows that 

[t ^ (Xt)*] , [t ^ {Xt),] G £{B\JC^+") (4.10) 

where we denoted the pull-back and push- forward with Xt by (Xt)* and (^t)*, re- 
spectively. Moreover, these operators are uniformly bounded in norm in a small time 
interval. 

Theorem 4.1. Assume that v G BUC^+" be given with v{0) = g e BUC^+". 
Then, forT > small enough, there exists a unique solution s of (|4.1|l - (|4.2|l on [0,T] 
with 

s G C^+^ BUC^+" n & BUC^+" . 

The existence interval is independent of v in a neighborhood of g. Furthermore, if 
vi,V2 G C'' BUC'^^" are such that vi[0) = W2(0) — g, then 

\\si - S2llci+''BUC2+°nCf BUC3+° 

< CT -U2||c0BUC3+° +clki -W2l|cfBUC2 + ° • (4-11) 

Proof. The solution can be computed by solving the reduced system of charac- 
teristic equations H4.8|l - H4.9|l . which produces in particular the diffcomorphisms Xf. 



18 



P. GUIDOTTI 



A minimal interval of existence can be chosen independently of the initial conditions 
thanks to the assumptions on v. Next 

shows together with the discussion preceding the formulation of the theorem, in par- 
ticular H4.10|l . that 

The desired time regularity has to be backed out from equations (|4.5|l - H4.(j|l . The time 
derivative appearing in the right-hand-side of H4.6|l is obviously not welcome in view 
of the assumptions made on v. It is, however, possible to rewrite the term containing 
said time derivative while integrating the equation to read 



r{t,x) = y/1 + \p{t,x)\'^v{t, x)+ 

, ,vV,vp]{T,XAXr'i^)))dr^ f\^fl=v]{T,Xr{X-\x)))dT. (4.12) 

The solution s can then be obtained by one further integration from (|4.5|) . By H4.10|l 
and the regularity assumption on v it then follows that 

The use of this argument needs to be justified. This can be done by substituting v 
by a regularized version of it which is differentiable in time and for which the partial 
integration used to obtain (|4.12|) can safely be performed. Letting it converge back 
to the original function v produces a function satisfying the modified equation. Since 
the solution of the original problem can be constructed by solving only the reduced 
problem H4.8|l - (|4.9() no trouble is encountered in taking the limit. 
The additional estimate follows from the characteristic system (I4.3(l - H4.7|l . the regu- 
larity assumptions on vi, V2 and the fact that vi(0) — V2{0) = 0. □ 

5. Existence Result. The local existence result for regular solutions of Hl.l|l - 
1)1.5(1 is based on the analysis of H1.6|l - H1.10|l . The regularity theory for linear singular 
elliptic and parabolic problems described in the previous sections and the analysis for 
the Hamilton- Jacobi equation performed in the last section will be the main tools. It 
should be kept in mind that e = 1 without loss of generality. To better highlight the 
structure of the problem and the use of the linear regularity theory it is convenient 
to rewrite (|1.6I) - (|1.8() using (|1.9|l in the following form 

u~A{t)u^A{s)u in 5 (5.1) 
u = g on To (5-2) 

1 

T 



-dyU = H{s,u) on Fi (5-3) 



where A{t) = A.j. + -p^rj-pjdy + ^dy and the operator A{s) is defined by 

t'^g'^{x)^ ^ ' "^s " s 



sAs-2|Vs|2 

y 2 ^y'^ 
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and iJ(s, u) is given by 

si si 

H(S,U) = : rr(Vs Vw) u(l + u) , 

^ ' t l + |VsP I ' t ^l + |Vs|2 ^ 

respectively. 

Remark 5.1. It should be observed that the operator family A{t) does not coincide 
with the one considered in the sections devoted to elliptic and parabolic regularity 
theory and denoted by A{t). It is, however, easy to check that is enjoys the same 
"singular" elliptic regularity theory as a perturbation argument shows 

A{t)^A{t)[idE,+A-\t)^dy] = [idE, +^dyA-\t)]A{t) 

Associated with the new family, there will be boundary solution operators Roit) and 
Rnit) with corresponding mapping properties. These new boundary operators are 
actually given by 

Roit) = Roit) + A{t)-^[^dyRD{t)] and RN{t) = RN{t) + A{t)-^['^dyRN{t)] , 

where Rd and Rn are the original boundary operators. From now on the old notation 
will apply to the new operators. 

Assume now that s e §, that is, that s e 0^+^^ BUC^+" n & BUC^+" is chosen in a 
small ball about the function [t > tg(x)] (in the natural norm of §) and satisfying 

s(0)=0, s(0) = .g. (5.4) 

It can be verified that the operator A{s) G £(Ei,Eo) is small in the operator norm 
if the time interval length T > is chosen small (this the reason why the modified 
operator A{t) is introduced, see Remark I5.1|l . This structure will be useful in the 
proof of the following existence result. 

Theorem 5.1. Let s G S 6e given with the above properties (|5.4(l . g G BUC^^", 
and let T > be small. Then there exists a unique solution u G Ei of H5.1|I - (|5.3II 
which is the fixed-point of the operator 

$i(u) ^ V + RD{t)g + RN{t)H{s,u) 

where v solves 

V ~ A{t)v = A{s)u ~ ^A-i(i)[A, + ^-dy\ [i?D(i)ff + RN{t)H{s,u)\ (5.5) 
Moreover, for S\,S2 G §, the following estimates 

II71U1 - 7iM2||c'3buc3+^^ < c||si - S2IIS/3 , 

II71W1 - 7i"2||c''BUC2+/3 < cT^llsi -S2IIS/3 (5.6) 

hold. 

Proof. As already pointed, the existence proof relies on maximal regularity results 
obtained in the previous sections which are at the heart of the matter and make it 
now possible to use a simple Banach fixed-point argument. The fixed-point u will be 
looked for in the set {u G M^-^{g, r) \ u{Q) ~ g} for positive r to be fixed later. This set 
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is endowed with the Ei topology and is therefore complete. It is easy but tedious to 
show that V, R]y{t)H{s, u) are in Ei n Eo_o- The first regularity claim follows directly 
from Theorems 13.31 and 13.41 both for v and RN{t)H(s,u). The vanishing property 
for V follows from its regularity using the equations it satisfies. Its regularity readily 
implies that dyv{t = 0) = 0. Since v satisfies homogeneous boundary conditions 

«(y = 0) = = dyviy = 1) 

it then is v{t = 0) = 0. As for the boundary term R]\[{t)H{s, u) it follows from 

combined with 

tH{s,u) e C|?BUC^+", 

or, alternatively, also from the equations. Finally the equations also imply that 

d^RD{0)g = , RD{0)9{y - 0) - .g , dyRD{0)g{y = 1) = . 

It clearly follows that Ro{Q)g = g. Summarizing it is obtained that $i(u)(0) = g. 
Next a norm estimate for — g is needed. The term RD{t)g — g can be estimated 

by some fixed constant ri. As for the other boundary term one has that 

\\^^[tH{s,u)]\k, < cr{l + r)T^\\s\k 

and can therefore be made arbitrarily small by reducing the interval length. The above 
estimate is possible since the map [t ^ tH{s, u)] vanishes at the origin in view of the 
properties of s. Next v has to be estimated. Most terms can simply be estimated in 
a way that they produce a factor and can be made arbitrarily small by interval 
length reduction. This is due to the fact that 



\v\\e, < \\{dt~A{t)) '\\cm)\\F{s,u)\\e, 



where F(s,u) summarizes the terms in the right-hand-side of (|5.5|l . and the specifics 
of F. Consider for instance the first term A{s)u. Since A{s) vanishes at the origin, it 
can be seen that 

\\A{s)\\ciE,,Eo)<cTPms 

which implies the claimed estimate. The estimate for the second term follows from 
the fact that 

and the general mapping properties of the operators involved. All remaining terms 
but one can be handled similarly. All terms in ^[tH{s, u)] containing either s or Vs 
are going to lead to more terms vanishing at the origin. The term 

^[ , ' u{l^u)] 
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behaves slightly difTerently and needs separate consideration since su{l+u) / a/I + |Vsp 
does not vanish at the origin. Rewriting this term as 

and using that u(0) = g up to vanishing terms of type , it follows that 

II ^ [ , u{l + u)] IImo <r, + cr- T%h 

where the constant r2 is given by 

r2 = ||(9,-^(i))^^:^[.g^(l+.9)]|k. 

It follows that $1 is self-map if r is chosen to satisfy r > ri + r2 and the time interval 
length is sufficiently small. It is important to observe that all estimates have constants 
which are independent of s in the chosen neighborhood. 

Next it needs to be shown that <&! is a contraction. The calculations are very similar 
to the above. The following terms 

- $1(U2) ^Vi-V2+ RN{t) [H{s, Ml) - H{S, 1*2)] 

need to be estimated. The difference vi — V2 satisfies 

Vi-V2 = {dt - A{t)y' [Fis.ui) - F{s,u2)] . 
Maximal regularity yields that 

\\vi-V2\\e, = \\{dt - A{t)y'[F{s,ui) - F{s,U2)]\\e^ < c\\F{s,ui) - F{s,U2)\\eo 

so that only the last term needs to be estimated. As for the other summand, an 
estimate for 

\\RN{t)[H{s,Ui) - H{s,U2)]\\E^ 

needs to be established. Using similar estimates as for the self-map property and 
observing that all the terms not vanishing at the origin drop out since they do not 
depend on the unknowns one arrives at 

||$l(s,Ul) - $l(s,lt2)||Ei < cr'^||wi-U2||Ei 

with a constant which is independent of s in the chosen neighborhood. Existence of 
a unique fixed-point follows by making the interval length short enough. 
Next observe that any solution actually satisfies 

7iw ^ -iiRD{t)g + 7i[« + RN{t)H{s, u)] 6 & BUC^+" (5.7) 

by Remark |2 . II since g G BUC'^+". This follows from 

v={dt- A{t))-\F{s,u) + li^^.(^).g2(l +g)] _ (a, - A{t))-'-^RN{t)g\l + .9) 
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and 

RNit)His, u) = RNit)[His, u) - 5^(1 + g)] + + g) , 

where the first terms on the right-hand-sides above belong to Ei^o and the second 
ones are not singular due to the regularity of 5, and therefore of +5), combined 
with Remark |2.1I The first continuous dependence result follows along the lines of 
the above estimates from the fact that the the solution to the initial boundary value 
problem depends smoothly (linearly) on the interior and boundary data A{s)u and 
i?(s, u). The only difference stems from the fact that a stronger norm is now estimated 
which leads to the absence of the factor T*^ in front of ||si — S2||s- This factor can 
be regained if a weaker norm is estimated in view of the "desingularizing" properties 
of [dt — A{t)]^^ , which leads to the stated continuous dependence estimate. The 
estimates are based on 

II71U1 -7iU2||c-3BUC3+'3 < ||7l"J'l(si,'"l) - 7l^'l(s2,U2)||c'5BUC3+'3 
< ||7l"I>l(si,Ul) - 7i'I>i(si,W2)||c'^BUC3+'' 

+ ||7l*l(si,W2) - 7l*l(s2,'"2)|lcf BUC3+'5 

The estimate for the first term of A{si) ~ .4(52) is typical. It gives 

II [4(1 + y'lvsip) - 4(1 + y'l v.2n] f k(E„Eo.o) 

< l|-(l + J^'lVsip) - -(1 + y2|Vs2|')||c^ BUCi + = C l|-7|-||£(Ei,Eo) 

Si S2 r 

which easily entails the claim since 




< c(pi - sallc-^ Buci+° + Pi - s2||c/Jbuc2+°) < c||si - S2II5 ■ 
One is eventually lead to 

||7l"l - 7l"2||c''BUC3+° < cT'^Wm - M2IIE1 + ||si - S2\\sfi 

from which the claim follows since the first term on the right can be absorbed in the 
left-hand-side. □ 

All pieces are now in place in order to show existence for the original system of 
equations 6(1 - ^.10(1 . By denoting with u = ^i{s) the solution of p. 6(1 - 1(1. 8(1 for a 
given s G 5, it follows that 

u e El and 7iu £ BUC^+" . 

Decomposition (15.7(1 shows that it is in a given neighborhood of g for any choice 
of s e 5 in a neighborhood of [t t-^ tg{x)]. Similarly let ^2{u) be the solution of 
((1.9(1 - 1(1.10(1 constructed in section|3| In this case, if u e Ei with 71U e C''BUC^+", 
^2{u) £ S will be in a neighborhood of tg{x) for u in a neighborhood of g. It follows 
that a solution to the full problem can be found by producing a fixed-point s for the 
map $ (f>2 o *l'i and defining u :— ^i{s). Summarizing 
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Theorem 5.2. For any given g G BUG such that g{x) > go > , x <E W\ 
there exists a unique local solution (w, s) of the free boundary problem (|1.6|) - (|1.1U|I 
such that 

u^ii^d'^^d^u G C^BUC^+"C /or < |a| + /c < 2 

7iM e BUC^+" and 

s e C^+'^ BUC^+" n C'^ BUC^+" . 

Furthermore 

u = RD{t)g + RNit)H{s, u)+v 

where v is a solution of (|5.5|) and the different terms have different asymptotic behav- 
ior in the origin. They behave respectively like t^ for fc = 0, 1,2 modulo the singular 
behavior built-in in the space chosen. 

Proof. Making the time interval as small as needed and choosing r > large 
enough it can be seen that $2) is a self-map on the complete set 

{u e 1ei (0, r) I u{Q) = g} x {s e Mgt, go/2) \ s(0) - , s(0) = g} . 

where the additional requirements u{0) = g , s(0) = , s(0) = g are included in the 
definition of the balls. Combining the estimates of Theorems 14.11 and 15. II it follows 
that 

||<f(si)-$(s2)||s. <CT'^||S1-S2||s^ 

and the contraction principle can be applied to obtain existence. □ 

Remark 5.2. It should be pointed out that this is the first well-posedness result for 
this class of singular free boundary problems in more than one space dimension for the 
full evolutionary problem. A companion result for the quasi- stationary approximation 
(e = 0) has previously been obtained in /?/. 



